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. ABSTRACT 

^ Context. X-ray observations revelled a group of radio-silent isolated neutron stars (INSs) at the centre of young supernova 
$^ remnants (SNRs), dubbed central compact objects or CCOs, with properties different from those of classical rotation-powered 
pulsars. In at least three cases, evidence points towards CCOs being low-magnetized INSs, born with slow rotation periods, and 
, possibly accreting from a debris disc of material formed out of the supernova event. Understanding the origin of the diversity 
of the CCOs can shed light on supernova explosion and neutron star formation models. Optical/infrared (IR) observations are 
crucial to test different CCO interpretations. 
^ : Aims. The aim of our work is to perform a deep optical investigation of the CCO RX J0822. 0—4300 in the Puppis A SNR, one 

■ of the most poorly understood in the CCO family. 

J — , Methods. By using as a reference the Chandra X-ray coordinates of RX J0822. 0—4300 we performed deep optical observations 
f — . in the B, V and I bands with the Very Large Telescope ( VLT). 

' Results. We found no candidate optical counterpart within 3a of the computed Chandra X-ray position down to 5a limits of 

, B~ 27.2, V ~ 26.9, and I ~ 25.6, the deepest obtained in the optical band for this source. 

■ Conclusions. These limits confirm the non-detection of a companion brighter than an M5 dwarf. At the same time, they do not 
constrain optical emission from the neutron star surface, while emission from the magnetosphere would require a spectral break 

, in the optical/IR. 
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1. Introduction the case, e.g. for the Crab pulsar, but have strong ther- 
mal components. Furthermore, CCOs show no evidence of 

Observations performed with the Etnstein satellite lead to ^^^^^^^ ^-^^^ ^^^^1^^ (p^^^^^ X-ray pulsations have been 

the discovery of puzzling X-ray sources in young 2-10 detected for two CCOs, IE 1207-5209 (P=424 ms; 

kyears) supernova remnants (SNRs), i.e. IE 1613-5055 in ^^vhn et al. 2000) and CXOU J185238.6+004020 in Kes 

RCW 103 (Tuohy & Garmire 1980) and IE 1207-5209 in ^p^^g^ ^^^^^^^^ 2005). Interestingly, in both 

PKS 1209-52 Helfand & Becker 1984 . Their location, ^ , r ■ 

^ ' ' cases the measured upper limits on the period deriva- 

very close to the geometrical centre of the remnant, siig- ^^^^^ ^^^^^^^ ^pj^ ^^^^ ^^^^ ^^^^^ ^ ^^^^^^ ^ ^^3 

gested that they were young neutron stars formed by the ^^^^^^ ^^^^ g^j^ ^g^^ ^ ^ ^jp^l^ magnetic fields 
supernova explosions, despite no pulsations being detected ^ ^ 5 ^ ^ rotational energy losses E < lO^^-ss 
at radio and X-ray wavelengths. About ten of these X-ray ^^-2 ^-i (Cotthelf & Halpern 2007; Halpern et al. 
sources, dubbed central compact objects or CCOs (Pavlov ^^^^^^ possibihty is that these CCOs were born spin- 
et al. 2002), have been discovered to date. The CCO X-ray ^.^^ ^^^^^ ^^^.^^ p^^^^^^ p^^^^^ ^.^^ ^^^^ 
properties markedly differ from those of young rotation- ^^^-^ g^j^^^ enabling polar cap accretion from a debris 

powered INSs (see De Luca 2008 for a recent review), j-r i-tj-i, j.- 

^ , ^ . ! disc formed out 01 the supernova explosion. In the opti- 

Their X-ray spectra are not purely magnetospheric, as is w r j n-n\ j u j.- u r j 

r' .7 to i 7 cal/mtrared (IK), deep observations have been performed 

~ — ; — ~ ^ , , , only for a handful of objects, mostly resulting in non- 
Sena offprint requests to: rm2Qmssl.ucl.ac.uk i , , ■ i ■ i 1 1 ^ , , 1 r 
J, , , . „ , ^^r^ ^ , , detections, which do not rule out the presence of a very 
Based on observations collected at ESO, Paranal, under 

Programme 78.D-0706(A) 
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low mass companion star (M type or later), or of a debris 
disc (see De Luca et al. 2008, and references therein). 

One of the most interesting CCOs is that in the 
Puppis A SNR (3700 years). The X-ray source (RX 
J0822. 0-4300), only hinted at in Einstein data (Petre et 
al. 1982), was then clearly detected by ROSAT (Petre et 
al. 1996). The source X-ray spectrum is predominantly 
thermal (Hui & Becker 2006) and well fitted by a two- 
component blackbody model. The corresponding hydro- 
gen column density {Nh ^ 4.54 x 10^^ cm~^) is consistent 
with the 2.2 kpc distance derived from neutral hydrogen 
observations of the SNR (Reynoso et al. 2003), which yield 
a 0.5-10 keV X-ray luminosity of Lx ~ 2 — 5 x 10^^ erg 
s-i for RX J0822.0-4300. Recently, a period at - 112 
ms has been reported from XMM-Newton observations by 
Gotthelf & Halpern (2009). The first optical observations 
of RX J0822.0-4300 yielded upper hmits of B~ 25 and 
R~ 23.6 on the source brightness (Petre et al. 1996) while 
new optical/IR observations (Wang et al. 2007) did not 
reveal any candidate counterpart down to 26.5 and 
Ks - 20.1 

Here we present deeper optical observations of RX 
J0822. 0—4300, performed with the Very Large Telescope 
{VLT} of the European Southern Observatory (ESO). 
Observations and data reduction are described in Sect. 
2, while the results are presented and discussed in Sect. 3 
and 4, respectively. 

2. Observations 

2.1. Observations and data reduction 

Optical imaging observations of RX J0822. 0—4300 were 
performed in service mode with the VLT Kueyen tele- 
scope at the ESO Paranal observatory in November 2006 
and March 2007 (see Table 1). The latter observations 
were performed right after the recoating of the telescope 
main mirror on March 2nd 2007. We used the FOcal 
Reducer/low dispersion Spectrograph {FORSl), a multi- 
mode camera for imaging and long-slit/multi-object spec- 
troscopy. At the time of our observations the camera 
was still equipped with the original single chip, four port 
2048 X 2084 Tektronik CCD detector. We selected the stan- 
dard resolution mode, with a 0'.'2 pixel size, and a field of 
view of 6'8 x 6'8. Sequences of 560 s exposures were ob- 
tained through the Bessel B, V and I filters, for total in- 
tegration times of 14000, 2800, and 11200 s, respectively. 
In order to avoid light contamination at the target posi- 
tion, bright stars in the field of view were masked using 
the FORSl occulting bars. 

Observations were carried out with airmass better than 
~ 1.3, and ~ 0'.'8 seeing. Sky conditions were photo- 
metric (or close to) for all nights, with the exception of 
that of November 26th 2006. For this reason, and due to 
the limited amount of exposure time, in the following we 
do not make use of this observation. Bias, twilight flat- 
fields, and observations of standard stars from the Landolt 
fields (Landolt 1992) were obtained as part of the FORSl 



Table 1. Log of the VLT/ FORSl BVI observations of the 
Puppis A CCO (RX J0822. 0-4300). Columns report the 
observing date, the filter name, the total integration time, 
the average seeing and airmass in each data set. 



yyyy-mm-dd 


Filter 


T(s) 


Seeing 


Airmass 


2006-11-26 


I 


2800 


0.88 


1.09 


2007-03-18 


B 


2800 


0.98 


1.06 


2007-03-20 


B 


5600 


0.84 


1.07 




V 


2800 


0.71 


1.14 


2007-03-23 


I 


8400 


0.71 


1.22 


2007-03-24 


B 


2800 


1.22 


1.05 


2007-03-25 


B 


2800 


0.89 


1.05 



science calibration plan. For the V-band observations no 
standard star observations were acquired on the night of 
March 21st. For photometric calibration purposes we thus 
used the closest-in-tinie available standard star observa- 
tions, taken on March 10th. 

Data were reduced (bias subtraction, fiat -field correc- 
tion, pre- and overscan trimming) using the ESO FORSl 
data;eductionpipelinfl We applied the photometric cal- 
ibration using the night zero points available in the in- 
strument data quality control databaseH. Airmass correc- 
tion was applied using the extinction coefficients for the 
Paranal Observatory. For the March 21st V-band data set 
the zero point trending plots ensure that the used closest- 
in-time zero point of March 10th, i.e. measured after the 
telescope recoating, does not introduce systematic eflFects 
in our photometric calibration. We converted the database 
zero points, computed in units of electrons/s, in units of 
ADU/s by applying the corresponding electrons-to-ADU 
conversion factors. For each dataset, we finally used the 
tasks drizzle and imcombine to align and co-add 
all exposures taken with the same filter, as well as to mask 
bad pixels and filter cosmic ray hits. 

2.2. Astrometry calibration 

The astrometric calibration was computed on the FORSl 
B-band image, taken as a reference, using the position and 
coordinates of 50 well-suited stars (e.g. not saturated, or 
too faint, or too close to the CCD edges, and evenly dis- 
tributed in the field of view) selected from the GSC-2 ver- 
sion 2.3 (Lasker et al. 2008). The pixel coordinates of the 
GSC-2 stars were measured by fitting their intensity pro- 
files with a Gaussian function using the dedicated tool of 
the Graphical Astronomy and Image Analysis {GALA) in- 
terfacqj- The transformation between the detector and sky 
coordinates was then computed using the Starlink code 
which is based on higher order polynomials and 
accounts for unmodelcd CCD distorsions. The rms of our 

^ http:/ /www.eso.org/observing/dfo/quality/FORSl/pipeline 

^ http:/ /www.eso.org/observing/dfo/quality/ 

^ star-www.dur.ac.uk/~pdraper/gaia/gaia.html 

^ http:/ /star-www. rl.ac.uk/Software/software. htm 
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Fig. 1. 20" X 20" cutout of the co added FORSl B-band 
image (14000 s) of the RX J0822.0-4300 field. The CCO 
position is marked by the circle (0'.'7 radius, 3a confidence 
level). The two objects detected closest to the Chandra 
position, also detected in the optical images of Wang et 
al. (2007), are labelled. 

astromctric solution is Sr w 0"15, accounting for the rms 
of the fit in the right ascension and declination compo- 
nents, to which we added in quadrature the uncertainty in 
the registration of the FORSl image on the GSC-2 refer- 
ence frame. Following Lattanzi et al. (1997), we estimated 
it as \/3 X (tgsc/VNs, where the \/3 term accounts for 
the free parameters (x-scale, y-scalc, and rotation angle) 
in the astrometric fit, acsc = 0"3 is the mean positional 
error of the GSC-2 coordinates (Lasker et al. 2008) and 
Ns is the number of stars used for the astrometric cali- 
bration. The uncertainty on the reference star centroids is 
below O'.'Ol and was neglected. We also added in quadra- 
ture the 0"15 uncertainty (Lasker et al. 2008) on the link of 
the GSC-2 to the International Celestial Reference Frame 
(ICRF). Thus, the overall accuracy of the FORSl astrom- 
etry is 0'.'22 {la). 

3. Results 

As a reference for the RX J0822. 0—4300 position we 
used the Chandra coordinates reported in Winkler & 
Petre (2007), i.e. a(J2000) = 08''21™57.355", (5(J2000) = 
-43°00'ir.'17 (epoch 2005.32; MJD= 53485), with a la 
uncertainty of iO'.'Ol per each coordinate. To register the 
RX J0822.0-4300 position on our FORSl B-band image 
(epoch 2007.23; MJD=54184) we then accounted for its 
Chandra proper motion of 165 ± 25 mas yr~^ (Winkler 
& Petre 2007), obtaining a(J2000) = 08^21™57.309^ 
5(J2000) = -43°00'17'.'28. We finally estimated the over- 
all uncertainty on the RX J0822. 0—4300 position at the 



epoch of our FORSl observations by adding in quadrature 
the positional error on the Chandra coordinates at the ref- 
erence epoch (0'.'015), the coordinate uncertainty due to 
the proper motion extrapolation (~ 0'.'048), and the error 
of our astrometric calibration (^ 0"22). This yields to an 
overall uncertainty of 0"23 (la). We note that both the 
Chandra coordinates and the GSC-2 ones, used for the as- 
trometric calibration of our FORSl images, arc all linked 
to the ICRF so that we exclude any systematic offset in 
our astrometry. 

Figure 1 shows the computed Chandra position of RX 
J0822.0-4300 overlaid on the co-added FORSl B-band 
image. We note that our error circle is offset by ~ 1'.'3 from 
that shown in Figure 1 of Wang et al. (2007). This is due 
to the X-ray source proper motion between the epoch of 
our VLT observation (March 24th 2007) and the Chandra 
one (January 1st 2000) from which they derived the ref- 
erence coordinates for their Magellan observation (April 
5th 2003). No candidate counterpart is detected within 
a conservative ~ 3ct error circle. We have computed the 
flux of the field objects through customized aperture pho- 
tometry. The closest objects, marked as 1 (B=25.6) and 2 
(B=25.8) in Figure 1, are detected at 2'.'6 and 4'.'7 south- 
east of the CCO position, respectively. Their measured 
offests are well beyond any possible residual uncertainty 
in our astrometry, which rules out any association with the 
CCO. No candidate counterpart is detected either in the 
V-band image or in the I-band one. We thus conclude that 
the optical comiterpart of the Puppis A CCO is unidenti- 
fled. We estimated the detection limits of our images by 
extrapolating the fluxes of the faintest objects detected 
close to the RX J0822. 0-4300 position. To this aim, we 
run the object detection using the SExtractor program 
(Bertin & Arnouts 1996) and we computed their magni- 
tudes through customized aperture photometry. We thus 
derived 5(7 upper limits on the optical brightness of the 
RX J0822.0-4300 counterpart of B~ 27.2, V ~ 26.9, and 
I ~ 25.6. 

4. Discussion and conclusions 

Our BVI flux limits conflrm the non-detection of a com- 
panion star earlier than an M5 dwarf, as suggested by 
the previous optical and near/mid- IR flux upper limits of 
Wang et al. (2007). Similar conclusions were recently re- 
ported for other CCOs (see De Luca 2008 and references 
therein). The available upper limits can also be used to 
constrain the optical/IR spectrum from the neutron star 
or from a possible fallback disc. In Figure 2 we plotted the 
extinction-corrected near/mid-IR and optical spectral flux 
upper limits together with the extrapolations of the best 
flts to the X-ray spectrum computed using double black- 
body models with the hydrogen column density Nu left as 
a free parameter (Hui & Becker 2006). We computed the 
interstellar extinction correction in the JHKg and BVRI 
bands from the Nh derived from the X-ray spectral fits 
applying the relation of Predehl & Schmitt (1995) and the 
extinction coefiicients of Fitzpatrick (1999). For the mid- 
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Fig. 2. Extrapolation in the optical/IR domain of the X- 
ray spectrum of RX J0822. 0—4300. Both curves corre- 
spond to a double blackbody best-fit spectrum (Hui & 
Becker 2006a) with (i) Ti = 2.6i x 10^ K,r2 = 5-04 x 10^ 
K, and emitting radii R\ ~ 3.29 km, i?2 ~ 0.75 km (solid 
line) and (ii) Ti = 1.87x 10^ K, T2 = 4.58x10^ K, i?i = 10 
km, i?2 = 1.09 km (dashed line). The extinction-corrected 
optical/IR flux upper hmits (Wang et al. 2007 and present 
work) are plotted in black and red, respectively. 



IR bands the interstellar extinction was computed from 
the Ay using the relations of Chiar & Tielens (2006). For 
simplicity, we plotted only the extinction-corrected spec- 
tral fluxes corresponding to an N}i — 4.54 x 10^^ cm^^ 
derived from the first double blackbody spectral fit. As 
seen, the X-ray spectra extrapolations are at least two or- 
ders of magnitude below the optical/IR upper limits, even 
in the more favorable case of thermal emission from the 
bulk of the neutron star surface (dashed line in Figure 
2). On the other hand, for a blackbody plus power-law 
model (Ti = 8.57 x 10^ K; F = 4.29) the X-ray spectrum 
extrapolation falls about four orders or magnitudes above 
the optical/IR upper limits. Thus, we can not constrain 
optical/IR emission from the neutron star surface, while 
emission from the magnetosphere would obviously require 
a break in the source optical-to-X-ray power-law spec- 
trum. We note that the spectral break implied by the bro- 
ken power-law model fit to the X-ray spectrum (Fi = 3.61; 
F2 = 5.29) would also over-predict the optical emission of 
RX J0822. 0—4300 by an amount similar to the blackbody 
plus power-law model. However, detection of optical mag- 
netospheric emission would require an energy output and 
an optical emission efficiency comparable to those of the 
Crab pulsar. 

The presence of a fallback disc can not be ruled out by 
our data. For the assumed 0.5-10 keV X-ray luminosity of 
Lx ~ 2 — 5 X 10'^'^ erg s~^, the disc emission would be dom- 
inated by the reprocessed X-ray flux (Perna et al. 2000). 
The X-ray luminosity of RX J0822.0-4300 is within a 
factor of 2 from that of the CCO IWGA J1713-3949 in 
G347.3-0.5, i.e. Lx 2.8 x 10^3 erg s^ at a distance of 1.3 



kpc (see Mignani et al. 2008 and references therein). This 
would imply, for the same set of disc model parameters, 
i.e. the disc inner and outer radii, the disc accretion rate, 
and the viewing geometry (Perna et al. 2000), a compara- 
ble disc IR luminosity. Thus, the optical/IR spectral flux 
upper limits of RX J0822. 0—4300 are still consistent with 
the presence of, e.g. an undetected disc with outer radius 
smaller than k, IRq. 

The present upper limits on the RX J0822.0-4300 
optical emission can only be marginally improved with 
the VLT, or with other 8m-class telescopes. On the other 
hand, much deeper IR observations of RX J0822. 0—4300 
are still feasible with the VLT and could help to further 
constrain the presence of a fallback disc. 
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